Comment on "Density perturbations in the ekpyrotic scenario 
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In the paper "Density perturbations in the ekpyrotic scenario", it is argued that the expected 
spectrum of primordial perturbations should be scale invariant in this scenario. Here we show that, 
contrary to what is claimed in that paper, the expected spectrum depends on an arbitrary choice of 
matching variable. As no underlying (microphysical) principle exists at the present time that could 
lift the arbitrariness, we conclude that the ekpyrotic scenario is not yet a predictive model. 
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I. INTRODUCTION 

Recently, a model called the ekpyrotic scenario m was 
proposed as a possible competitor to the inflationary 
paradigm || on the basis that it could solve the problems 
of the hot big bang model and also produce an almost 
scale invariant spectrum of primordial scalar perturba- 
tions. Although the foundations of the model are still 
under debate S, we only shall comment on the contro- 
versial |Q, H [| claim, crucial for the ekpyrotic scenario, 
that the power spectrum is scale-invariant 0. 

In Ref. JtJ, the complicated five-dimensional evo- 
lution of the background is modeled by an effective 
four-dimensional spacetime which experiences a singu- 
lar bounce at n = (77 being the conformal time) as 
the effective scale factor vanishes. The calculation of the 
power spectrum is a controversial issue as different au- 
thors would use different matching conditions || |[ [u]] 
at the singular bounce. There, Bardeen's gravitational 
potential jflj <!> diverges. The authors of Ref. H sug- 
gested to work with the comoving density contrast e m 
which is regular at 7/ = and has two linearly indepen- 
dent modes D and E, such that e m = e D(rj) + e2E{rj). 
At the singular point the equality of the coefficients be- 
fore and after the bounce is assumed B: £o(0 + ) = eo(0~) 
and £2(0 + ) = e 2 (0~). Although this rule is given without 
justification (recall that the standard junction conditions 
follow from the Einstein equations) , we shall refer to such 
equalities as matching conditions hereafter. 

In Ref. p"2| ] we criticized the matching conditions sug- 
gested in Ref. J7j. We showed that they lead to an am- 
biguity, and thus an arbitrary choice needs to be made 
that present-day physics cannot make. Ref. was sub- 



sequently modified, see Ref. ]l3| , to address some critical 
comments that have been made on the ekpyrotic model, 
one of them being the above mentioned ambiguity ex- 
hibited in Ref. iQ. In the new version jl3|], more de- 
tails about the matching conditions are given and one 
claim has been added (claim 1): 11 One situation of spe- 
cial interest is ... where the potential is irrelevant at 
<f> — > —00 and there is no radiation in the incoming state. 
In this case, £2(0") = 0. In this case, we would obtain 
the same final result from any matching rule which set 
e 2(0 + ) = Ae2(0~), with any constant A." It concerns a 
particular case that we had not considered and for which 
our result supposedly does no t ap ply. It should be noted 
however that the rest of Ref. [[13| inconsistently relies on 
the general case instead of this particular one, as it was 
the case in the preceding version 0. Claim 1 clearly 
utilizes, in an essential manner, the ideas developed in 
Ref. and represents a tentative a priori response to 
the points Ref. |i2| raised. 

References ft |13| also state (claim 2): "The prescrip- 
tion is invariant under redefining the independent solu- 
tions, e.g. by adding an arbitrary amount of the solution 
E{rj) to D(n). Matching any other non- singular pertur- 
bation variable, defined to be an arbitrary linear combina- 
tion ofe m and e' m with coefficients which are non-singular 
background variables (defined to possess power series ex- 
pansions in t, as above) will, with the same prescription 
of matching the amplitudes of both linearly independent 
solutions, also yield precisely the same resulf . This is 
an important point in the ekpyrotic scenario because, if 
true, it would support the use of the variable e m as a tool 
to define matching conditions. 

Below we show that both quoted claims are incorrect. 



II. IS THE EKPYROTIC SPECTRUM UNIQUE? 
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Let us start with claim 1, i.e. the question concern- 
ing the invariance of the spectrum under rescaling of the 
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parameter £2 and let us first define our notation. The au- 
thors of Ref. jl3) consider the equation of state close to 
the singular bounce to be w = p/p = ujQ+ujin+^V 2 ^' ' '■ 
Before the bounce, the kinetic energy of the scalar field 
dominates and the dynamics can be approximated by 
that of a free scalar field for which ljq = 1 and u>\ — 
L02 = (note that one can show that the values of Wj, 
Vi > 2, are irrelevant for this discussion). This should 
apply for 77 < close to the singularity (|^| <C 1) and 
corresponds precisely to the situation that is assumed in 
claim 1 and to the case where the equations of Ref. |l2] | 
cannot be applied (as mentioned in that article), hence 
the argument of Ref. Below, we complete the ar- 

guments of Ref. jl2| to include this particular case and 
show that the conclusions of Ref. j[2| remain unchanged. 

At leading order, the behavior of the scale factor is a 
simple power-law, 0(77) = io(— r/) 1//2 /(2\/2), in agreement 
with the notation of Eq. (24) in Ref. iQ. The comoving 
Hubble rate and the sound speed are Ti = a' /a = l/(2n) 
and c 2 = 1. The solution for the Bardeen potential flO, 
[l4| in the long wavelength limit reads 



<I> 



" dr 



(1) 



with 6 = l/(a\fl + u>) and the integral sign stands for 
the primitive of the integration kernel. The coefficients 
B\ and B2 are functions of the comoving wavenumber fc, 
to be determined by means of a matching with the initial 
vacuum condition ffl, [l^] . In the case at hand, 



f 2 n 2 
*0'/ 



Bi + ~B 2 . 



(2) 



Note that the equation of motion for the Bardeen poten- 
tial for a free scalar field reads <&" + (3/77)$' + fc 2 $ = 0, 
which can be solved exactly in terms of Bessel functions, 



^FMii B2Jl( ^ )+ 2^ liVl( - fc??) J' (3) 

whose expansion in powers of 77 permits to recover the 
special solution (43-44) of Ref. (?]] as well as Eq. (||) . 

Before the bounce the density contrast on comoving 
hypersurfaces, e m = — 2fc 2 $/(3H 2 ), is given by 

8k 2 



k 2 V 2 B 2 . 



(4) 



Comparison of Eq. (g) with Eqs. (43-44) of Ref. @ yields 



e<(k) = ^B<, 



: (k) 



-k 2 B<, 



(5) 



showing that the two "modes" decouple. It is interesting 
to compare with the case uif ^ 0, for which one has |l^] 
(here we assume = for the sake of simplicity; this 
does not change in any way the conclusion) 

8/c 2 < 128 In 2 k 2 
-1T B i 



e 2 



9^ 
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(6) 



,<2 



-Bf - (l + 21n2)fc 2 J B 2 <. (7) 



Let us note at this point that, with the choice we made 
of the normalization of the integral in Eq. (|l|), the limit 
lui — > is singular, as can be seen from Eqs. (||) - (Q) 
and as already mentioned in Ref. p^| . Another choice 
is possible as one is free to add an arbitrary constant 
amount of B 2 into B±, simply by making the integral 
a definite one. In particular, one can arrange that the 
constant factor (In 2) in the above equations be made to 
vanish from the outset, and it is easy to convince oneself 
that such a choice does not modify anything if the correct 
equations are used. 

In the ekpyrotic scenario, and in the long wavelength 
limit, the coefficient Bf oc fc~ 3 / 2 is dominant with re- 
spect to B2 oc fc -1 / 2 . This means that the relation 
£2 — 9w 2 eo/64 holds true if uj\ ^ 0, whereas it does not if 
w\ = 0. In the latter case, one has £2 = — ■ffiB2Co/(8i?i ) oc 
ke ^ 0, see Eq. (||). The authors of Ref. Q missed that 
the limit u>\ — > does in no way imply £2 = exactly, as 
implied by claim 1, in which this value is said to be used 
to perform the matching. Moreover, as demonstrated 
below, the fact that £ 2 becomes subdominant does not 
remove the ambiguity noticed in Ref. p2| . 

This is also related to the excessive claim in [[[3) that 
11 there is no long wavelength contribution to £ in the col- 
lapsing phase" . Inserting (^|) into the gauge-invariant 
variable C = (2/3)(7i~ 1 $ / + $)/(l + uj) + $ (see e.g. 
@, 00), gives ' 



C - \B2M-kn) + l^Noi-krj) 



which, in the limit 77 — > 0, yields 

c 



X -B2 
2 2 
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- 7J B 1 [ln(-fc77) + 7 E ], 



(8) 



(9) 







where 7 E is Euler constant coming from the expansion 
of the Bessel function. Since the second term is singular 
the first term is subdominant, but that does not imply 
that it vanishes altogether. 

After the singular bounce (77 > 0), 7^ and one 



recovers the relations |12 
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8k2 B> 



128 In 2 k' 



9lo> 2 



(10) 



B> (1 



21n2)fc 2 B 2 > - (11) 



The authors of Ref. |13j propose to glue the epoch after 
the bounce to the epoch before the bounce by imposing 
that £0 and £2 are the same before and after the singular- 
ity. Contrary to standard junction conditions, and until 
some more microphysics is specified, this recipe does not 
rest on any physical principle and one may wonder what 
is the rationale behind it. The reason why one is pushed 
to adopt a new rule in this case is linked to the fact that 
one is dealing with a perturbative approach around a 
singular background whose meaning is questionable [H. 
Moreover, we argued in Ref. jl2| that this prescription 
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is ambiguous since an arbitrary rescaling factor f(u>,k) 
shows up in the final result. Therefore, the scenario con- 
tains an arbitrary function in a physically measurable 
quantity: until this function can be calculated unambigu- 
ously by first principles, the model cannot be falsified. 

As in Ref. we rescale £2 by the completely ar- 

bitrary factor 1/f. This gives the constant part of the 
gravitational potential 



9lu 



>2 

1 R< 



(12) 



This is the main equation of Sec. II. It represents the 
equivalent of Eq. (54) of Ref. [|l2| for the case ujf = 0. 
On the other hand, the decaying mode amplitude reads 



= -(l + 21n2)B 



16 In 2 ll 

9lu> 2 



P 

* 2 F 



(13) 



given here for the sake of completeness. 

Note that if the expansion of the equation of state pa- 
rameter is done to higher order in r), all the relations of 
this section are only changed by inclusion of UJ2 through 
the replacement (3/8)w? — > = uj 2 + (3/8)u; 2 , in 

agreement with Ref. [i], |l3| , while none of the Wi for i > 2 
does contribute. 

Eq. (|l^) shows that the spectrum effectively acquires 
a scale invariant piece oc B^, together with an arbitrary 
piece oc B^ f > / f K ■ Indeed, as discussed in Ref. [ GjJ , 
the functions f < and f > can depend on the background 
constants (e.g., uS) and on the considered scale k in a way 
which is not yet given from first principles. The choice 
/ > = f < fixes the k— dependence of the spectrum, but, 
without physical justification, this remains an arbitrary 
choice, equivalent to assume a scale invariant spectrum 
from the outset. 



III. AN EXAMPLE 

Let us now consider claim 2 according to which the 
choice of e m is essentially unique. The ambiguity in- 
herent to the method proposed in Refs. Jf), [TJ| can be 
most clearly exhibited by means of an example suggested 
at the Euro-conference in Annecy in December 2001 by 
G. Veneziano [[[5]. This section relies completely upon 
his idea. 

In Ref. [|l3| it is argued that the density contrast e m is 
the quantity of interest because it is finite at 77 = and we 
have shown above how this quantity is used to propagate 
the spectrum through the singularity. However, e m is 
not the only finite, physically relevant, quantity. Based 
on the method developed in Ref. Q, we could equally 
well consider and use the conjugate momentum II to the 
conserved quantity £ ||] which obeys the Hamilton like 
equations, valid for isentropic perturbations, 



n = z 2 C, n' = -fcWc, 



(14) 



where z 2 = a 2 (l + uj)/c 2 . It is worth pointing out that 
since the quantity £ diverges at the singularity, it may be 
more useful to combine Eqs. ([l4|) into the single one 



n" 



2] nw 



k 2 c 2 n : 



0. 



(15) 



In the neighborhood of the singularity, c 2 = to = 1, 
and the solutions of Eq. (15) are — kr/Ji(— krf) and 
— krjNi(—kri), where J\ and N\ are Bessel function of 
the first and second kind respectively. These solutions 
are completely regular at 77 = 0. Note that Eqs. (|T^ ) 
and ( |l5| ) apply precisely in the two cases of interest here, 
namely that of purely hydrodynamical perturbations and 
that of a free scalar field, by replacing c 2 — *■ 1 || in this 
last case. The quantity n is related to e m as 



n = a z He n 



(16) 



Note also that the use of n may be argued to be more 
appropriate in view of the fact that in the regular bounce 
case |fl7| , even though admittedly a different case as 
the ekpyrotic model, the variable e m ends up being an 
odd function of ry, whereas n is even. In the short- 
duration limit, that would imply that e m experiences a 
jump whereas n is continuous. 

Since \a 2 Tt\ ~ ^g/16 is finite at the singularity 77 = 0, 
this relation clearly shows that the variable n, being a lin- 
ear combination of e m and e' m (with vanishing coefficient 
for the latter), satisfies all the requirements of claim 2. 
Accordingly, it seems that there is no convincing reason 
to use e m rather than n. In the vicinity of the bounce, 
n reads 



n 



16 



e - -wie ?7- 



(2 



8 



. (17) 



where s = — 1 before the bounce and s = +1 after the 
bounce. The proposal of Ref. Q] then would consist in 
assuming that the coefficients in front of the constant 
term on the one hand and in front of the rf term on the 
other hand are the same before and after the bounce. In 
the present context, this reduces to 



3wW> 



(18) 



since u/ 2 -^ = by definition. From these relations, it is 
easy to establish that 



Rk 2 

— R< f > - k 2 R< 

l 



3wW>k 2 < 

~1T~ 1 



(19) 



Finally, the computation of shows that the spec- 
trum no longer contains a scale invariant piece, 



B> = -B< cx k- 1 ' 2 . 



(20) 



but gives a spectral index equal to 3. It should be no- 
ticed that the same calculation goes through, with the 
same resulting spectral index if, in Eq. (Jlq) , one replaces 
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the factor a 2 TL by its magnitude \a 2 TL\ which is not only 
regular but also symmetric across the bounce. Therefore, 
the choice of the variable used to apply the rule proposed 
in Ref. [Q plays a crucial role. In the absence of an un- 
derlying reason to choose a variable rather than another 
at the present stage, II seems to be as relevant as e m , 
and we are led to the conclusion that claim 2 is definitely 
erroneous. Moreover, II is not the only well-defined vari- 
able that leads to such a conclusion: there is an infinite 
number of such possible combinations. 

IV. CONCLUSIONS 

To conclude, let us mention yet another possibility. In 
Ref. |H][| , it was suggested that the condition used in 
Ref. rajmay be cast in the form of a matching of the en- 
ergy perturbation as seen by an observer comoving with 
the fluid, and that such a matching is "at least as nat- 
ural as matching the energy in the longitudinal gauge" . 
Firstly, one should notice that no rigorous proof exists 
that the conditions of Ref. fllOf are equivalent to the pre- 
scription utilized in Ref. J7|. Secondly, we have shown 
in Ref. |l2j that, in the well known case of radiation to 
matter transition for which one can compare with the ex- 
act solution, this matching condition yields an incorrect 
result. A similar conclusion has been reached in Ref. || 
in the case of the reheating transition in which the equa- 
tion of state also jumps. However, it can be argued that a 
bouncing situation is in no way comparable to a jump in 
the equation of state and that the usual matching condi- 
tions may then not apply. This is indeed what was found 
m Ref. It 71 for the case of a non singular bounce. The 
case for the ekpyrotic scenario is still more involved, as it 
was also argued that the bounce must be singular if the 



brane collision is to produce radiation |18| . In any case, 
matching through an unspecified bounce is still an open 
question as there is no geometrical or physical argument 
imposing some particular choice, and no conclusion can 
be made about the power spectrum of perturbations in 
general, unless some non-singular bounce is specified, as 
it is the case of Ref. |l7j . 

Finally, we would like to insist on the fact that the 
overall gauge-invariant perturbation theory may be com- 
pletely meaningless when a singularity is reached, since 
some gauge transformations that are admissible in any 
other situation could turn out to be singular. In such a 
situation, a "gauge-invariant" variable will, at the singu- 
lar point, become "gauge-dependent", and its use rather 
than the use of any other variable become an arbitrary, 
physically meaningless, choice. In fact, the transforma- 
tion between the zero shear (conformal Newtonian) slic- 
ing and the comoving slicing is indeed singular at the 
bounce, which actually means that one or both slices are 
unphysical. Lyth Q has shown that there is no slicing 
on which the density contrast and the intrinsic curvature 
perturbation are finite simultaneously. In such a situ- 
ation, linear perturbation theory becomes meaningless. 
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